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Abstract The aim of the present work is to introduce a method based on 
Chebyshev polynomials for the numerical solution of a system of Cauchy type 
singular integral equations of the first kind on a finite segment. Moreover, an 
estimation error is computed for the approximate solution. Numerical results 
demonstrate effectiveness of the proposed method. 
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1 Introduction 

Let us consider a system of singular integral equations of the form 



dr + [ K(t,T)<P{T)dT = F{t),—1 < t < 1 (1) 


J-i 


where 


K{t,T) = [Kij{t,T)], i,j = 1,2,...,N 

Fit) = [hit), hit), , fNit)f 
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^{t) = [01 (i), 02(i), ■ ■ • , 0jv(t)]^ , 

A{t) = [a^jit)] , i,j = 1,2,...,7V, 

B{t) = [bij{t)] , i,j = 1,2,...,N . 

Here, (^, and {/i(i)}^i are given real-valued Holder functions 

and {0j(t)}jLi are the unknown functions. The matrices A and B are known 
such that S = A+B and D = A —B are nonsingular for all f G [—1,1]. In some 
familiar physical problems the entries of the matrices A and B are constants. 

The singular integral equations play important roles in physics and the¬ 
oretical mechanics, particularly in the areas of elasticity, aerodynamics, and 
unsteady aerofoil theory. They are highly effective in solving boundary value 
problems occurring in the theory of functions of a complex variable, potential 
theory, the theory of elasticity, and the theory of fluid mechanics. A general 
theory of the system of equations o has given in m- 

We study the system m in the case that A{t) = 0 and B{t) is a constant 
matrix. Therefore, the equation of system o takes the form 



Studies on this singular integral equation can be found in some literatures 
(See [TJ[2l[5l[6]). Chakrabarti and Berghe [2] proposed a method for solving 
Eq. ([2]) using polynomial approximation and collocation points have chosen 
to be the zeros of the Chebyshev polynomials of the first kind for all cases. 
Kashfi and Shahmorad [5] constructed another approximate solution of this 
equation using Chebyshev polynomials of the first and second kinds. Some 
other methods for solving this equation can be found in m- A convergence 
analysis of Galerkin and collocation methods for Eq. m has been given by 
Miel [TU] . 

A special type of Eq. m is the famous Cauchy singular integral equation 

f ^ ^ -l<t<l, (3) 

J-lT -t 

which has the following analytical solutions in four special cases based on 
boundedness of the unknown function 0 (t) at the endpoints of the interval 

[-1, 1] [2ll8lfT^ . 


Case 1. If the function 0(r) is unbounded at the endpoints r = ±1, then 


0(t) = 


Oq 


1 




\/l — 7r^\/l — J_i 

where oq is an arbitrary constant. 


t — T 


dt, —1 < T < 1 (4) 


Case 2. If the function 0(t) is bounded at the endpoints r = ±1, then 

r-l 


0(t) = - 


\/l — 


f{t) 


1-1 \/I - {t - t) 


dt , —I < T < I 


( 5 ) 
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a necessary and sufficient condition of existing this solution is: 



fit) 


dt = 0. 


( 6 ) 


Case 3. If the function is bounded at the endpoint r = — 1 and un¬ 
bounded at the endpoint r = 1, then 


^(r) 



1 - t fit) 
1 + tt — T 


dt, —1 < T < 1. 


( 7 ) 


Case 4. If the function (/)(t) is bounded at the endpoint t = 1 and unbounded 
at the endpoint t = — 1, then 


(^(r) 



1 +1 f{t) 
1 — tt — T 


dt, —1 < T < 1. 


( 8 ) 


More methods for solving Eq. ([3]) have given in [31[71 [T31[T3] . 

In the next section, we investigate approximate solutions for system m in 
the above four cases. 


2 Approximate solution 

To find approximate solutions for system o in cases 1,2,3,4, for v G {1, 2, 3,4} 
we set 


M 


to 


n i'^) - i^) ■= \ XI ’ j = 1,2,..., A 


and 


M 


K,j{t,T) ■=^-fijkit) Py,kiT) , = 1,2,..., A, 




where 


Tjix) -- 

= cos(j6»), 

iy=l, 

Ujix) = 

sin(0'+l)6») 
sin(e) ’ 

V = 2, 

Vjix) = 

cos(0-|-i)6») 
cos(f) ’ 

V = 3, 

Wj{x) = 

sin(0'+l)®) 

sin(f) ’ 

1^ = 4, 


(9) 


( 10 ) 
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are the Chebyshev polynomials of the first to fourth kinds and 

K{t) = 



in which x = cos{0). The roots of Chebyshev polynomials P^, m+i{x) are given 
by 




COS 


/ (2n-l)7r \ 
I 2(M+1) ) 


COS 


M+2 J ’ 


COS ( 

y 2M+3 J ’ 

_ cos ( 2 !^), 


1 / = 1 , 

i/ = 2, 

= 3, 

iy = 4. 


( 11 ) 


where n = 1,2,... , M + 1. These roots are used as the nods of Gauss- 
Chebyshev quadrature rules. 

Lemma. [S] The Chebyshev polynomials satisfy the orthogonality property 

'0, j, 

TT, i=j = 0, V=l, 

f, v = 1, (12) 


/■' K{t) 
/_1 


Pu,i{t)Pv,j{t)dt = < 


f Pu,j (x) 

1-1 Vl - T -t 


dr = TT < 


* = 7, 

jy = 

i=j, 

u = 

integral, we have 


v = l. 


v = 2. 

Wj{t), 

V = ?,, 


v = A. 


(13) 


Now we describe details of finding approximate solution in cases 1-4. 
Case 1. For 1 / = 1 the relations (©-(Uni) take the forms 


1 “ , 

(t)j{T) ~ ipijir) := XI > i = 1 


= 1,2,...,7V 


(14) 


1=0 

















Approximate solution of a system of singular integral equations 


5 


K,j{t,T) :='^ j^jkit)Tk{T), i,j = 1,2,...,N (15) 

k—0 

where Pji {j = 1,2,..., N, 1 = 0,1,..., M) are unknown coefficients and 
the symbol (y) ) denotes that the first term in the summation is halved. The 
functions 


7yfc(0 


2 

TT 


/■i K,,{t,r)n{T) 

y_i 


b j = 1,2,... ,fV, 




can be determined exactly, or may be approximated by using Gauss-Chebyshev 
quadrature rule, i.e 


2 M+l 

lijkit) ~ M 1 ^ ^l,s)Tk{xi^s) 1 

s=l 

where xi^s obtain from (El- 

Substituting from E1-(E1 into Eq. m and using Ei-ini for V = 1, 
gives the system 

N M N M 

EE bijf3jiUi-i{t) + ^EE' lijk{t)l3jk =, i = l,2,...,N. (16) 

j—1 1—1 j — 1 k—0 


If the given functions fi{t) and jijkit) are square integrable on [—1, 1] 
with respect to the weight function then they can be expanded as 

^M-l 


= l,2,...,iV, 


J lijkit) ^ Ui{t) , i,j = 1 

I y^it) - Efio" cu Ulit ), 1 = 1,2,...,N. 


where the coefhcients 

Gijki = f J^iVl- f^^ijkiOUiit) dt 


k = 0,l,... ,M 
(17) 




= 4^ fulfil \J T)Ui(t) Tkir) drdt 

i,j = 1,2,...,N, k = 0,l,...,M, 1 = 0,1,...,M-1 

= 4^l4iVT^mUiit)dt, i = l,2,...,N, 1 = 0,1,...,M-1, 


can be approximately determined from 

4 2 


- (M+ir ^r=l T,s=l 0 - X 2 ,r)Kij{x 2 ,r,Xi^s)Ul(x 2 ,r)Tk{xi,s) ■ 


- ^(J+1) l(l - xl^)Mx2,r)Ul{x2,r)- 


(18) 
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Using (113 in (HU) and linearly independence of {Ui{t)}^ yield 


M-l N 


M-1 N M 


M-l 


EE ^EEE' Gijkl^jkUl(t^ — 

1^0 j^l 1^0 j^l k^o 1^0 

which leads to the linear system 


N 

E 

i=i 


1 ^ , 

^ij/3j{i+l} + ^ GijklPjk 


k=0 


= ca, i = l,2,...,N, l = 0,l,...,M -1 

(19) 

for the unknown values Pjk {j = 1,2,..., N, k = 0,1,..., M). By taking 
arbitrary values for example for Pn, ... ,Pni, the remaining coefficients Pjk 
are uniquely found via the linear system (HU) which determine the elements of 
the vector function <P{t) via Eq. (HU). 


Case 2. We set = 2 in (I^- (ITU1) and substitute them in Eq. (HJ to get 

N M ^ N M 

2 


N M N M 

-EE bij — EE lijk{t)Pjk = -Mt) , i = l,2,...,N. (20) 


j=l 1=0 


j—1 k—0 


where we used the formulae (HU-iHSl)- Then, we expand the functions fi{t) 
and "fijkit) as 

(l^Mt)-J:f^o'G^JklTlit), i,j = l,2,...,N, k = 0,l,...,M 


\ 




i=l,2,...,N 


where the coefhcients are determined by 

Gzjki = i,j = 1,2,...,N, k,l = 0,l,... ,M, 

= \J (t, r)Ti{t) Ukir) drdt 


i = l,2,...,N, l = 0,l,...,M 


^ — 7r2 /-I dt , 

or approximated by 

Gijkl — (m+1}{M+2) X)r=l X)s=l ~ ^2,s)Tl{xi^r)Uk{x2,s) , 

^ ^(Af+1) h{xi,r)Tl{xi,r)- 

( 21 ) 

Using the last expansions in Eq. (EOl), returns the following linear system of 
equations 


^ Efcl Efc=0 GijklPjk — Cil, 


2 


i = 1,2,... ,7V, 1 = 0 


l:'^=l{-KPJ{l-l} + \Y.k=oG^JklPJk} = cu, f = 1,2,... ,7V, 7 = 1,2,..., 


( 22 ) 

M 
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for the unknown values Pji (j = 1, 2,..., iV, 1 = 0,1; • ■ •; Then the 
elements of the vector function <P{t) obtain from Eq. Q. 

Case 3,4. Proceeding by the same way as we did in cases 1,2, we get the 
linear systems 

N ( M 'j 

+ i = l,2,...,y 1 = 0,1,...,M(23) 

j=i I fe=o J 


and 


N ( M 

}> =c,i, i = l,2,... ,iV, Z = 0,1,...,M(24) 


fc =0 


respectively for v = S and v = 4, and then we determine the elements of 
corresponding vector <?(t) via ([5]). 


3 An estimation error and numerical results 

In this section, we describe an estimation error for the approximate solution. 
Let 

W{t) = [(pi{t) , ip2{t), , (pN{t)f 

be the vector of approximate solution of the system ([1]) and E{t) = (p{t) — <P{t) 
be the associated vector valued error function. Due to the approximation 
for A{t) = 0 the system ([T]) may be written as 

^ ^ dr + [ K{t,T)<l>{T)dT = F{t) + H{t), —1 < t < 1, (25) 


1-1 


T — t 


1-1 


where the perturbation term H{t) obtains from 

H{t)= [ —^ dr + [ K{t,T)<P{T)dT — F{t), —1 < t < 1. (26) 

J-i T-t 

Subtracting Eq. (I25p from Eq. o, yields a system of error equations as 

[ dr + [ K{t,T)F{T)dT = H{t), —1 < t < 1 (27) 

7-1 T-t 

which is solvable approximately like the system O- 

The following examples illustrate application of the method. 

Example 1. Let 


A{t) = 0, B{t) = 


, /i(t)=7r, f2{t)=2Trt 
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and find the solution of system o in case 1. 

By the above information the system m reduces to 

f f-i ~ + /y t(j)2{T)dT = TT, -1 < t < 1 

[ hi + hi T4‘i{T)dT + hi('^ + t)4>2{'r)dT = 27rt, —1 < t < 1 

( 28 ) 

since the matrices S = A + B = I 2 and D = A — B = —I 2 are nonsingular, 
therefore the solution of system (ESI) exists. The kernels Kij{t,T), K 2 j{t,T) 
(j = 1,2), and the functions fi{t), f 2 {t) are polynomials of degree at most 1, 
so we set 

</>! (t) := {Pjo To (r) + /3 ,1 T, (r) + /3,2 T 2 (r)} , j = 1,2 (29) 

and 

Kij{t,T) = 7yo(t)To(r) +7yi(t)Ti(T), i,j = 1,2 
ft{t) = CioUo{t) + c^lUl{t), 1 = 1,2 

where 

7iio(i) = 7iii(t) = Uo{t), 7i2o(i) = 7i2i(0 = 0) 

72io(0=0, 72ii(0 = C^o(i)> l 22 o{t) = 7221 ( 1 ) = C^o(i), 

Cio(l) = 1, Cii(l) = 0, C2o(l) = 0, C 2 l(l) = 1. 


Substituting these expansions into 
obtain 


and using dnii-disi), for 12 = 1 , we 

/3iit/o(l) + Pi 2 Ui{t) - ^loUiit) + ^nUoit) + \^ 2 oUi{t) = Uo{t) 

(30) 


^2lUo{t) + p22Ul{t) + \PllUo{t) + \p20Ul{t) + \p2lUo{t) = C/i(i) 

Then the linear independency of {17o(l), C/i(l)} implies 

|/3ii = 1 

— ■J/dlO + /3 i 2 + ^/d20 = 0 
\Pii + |/32i = 0 
\P20 + fd22 = 1- 

A nonunique solution of this system for the arbitrary values of /3io and /I 20 is 
given by 


PlO, Pll — I3 i 2 — iPlO — P 20 ) , 1320, /321 — — g, 


1322 — 1 — 2^20- 


For example, if /3io = (32o = 2, then j3i2 = P 22 = 0 and so we hnd from (1^ 

!t + 2 , , -It+ 2 


Mt) = 




02 (r) = 


VT^ 


(See figure [T] for the behavior of these solutions). 
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I &phi;_1 Plot-&phi;_2 Plot] 



Fig. 1 The plots of approximate solutions of Example 1 for M=2. 


Example 2. Solve the problem of Example 1 in case 3. 

In this case, we set 

•= \/t^ : j = 1, 2 (31) 

and 

Kij(t,T) = 7ijo(t)Vb(T) +%ji{t)Vi{T), i,j = 1,2, 

Mt) = CioWoit) + CiiWi{t), * = 1,2, 

where 

7iio(0 = Wo{t) - 7 iii(t) = 72io(0 = 72ii(0 = 722i(t) = 

li2oit) =-^Wo(t) + 7i2l(t) = 0, "f220(t) = ^Wi(t), 

Cio(t) = 1, Cii(t)=0, C2o(t)=-l, C 2 l(t) = 1. 


Substituting these expansions into (EHD and using (nm-disi) for 12 = 3, result 

f3wWo{t) + PiiWi{t) + l3ioWo{t) — \fiwWi{t), 

+\PiiWo{t) - i/320W'o(i) + \li 2 oWi{t) = Wo{t) 

< (32) 

/32obbo(0 + P2iWi{t) + \PwW(i{t) + i/3iiVlo(t) 

+\P 2 oWi{t) + \l32iWo{t) = -W'o(t) + Wi{t) 

and from the linear independency of {lEo(t), we get the algebraic 

system 

( 2/3io + ^/3ii — 4/320 = 1 

I —\l 3 w + /3ll + 2^20 = 0 

I + \Pll + P 20 + 5/321 = — 1 
i i/320 + /321 = 1, 


which has the unique solution 


/3io 


10 


/3ii 



22 


/321 


20 

y ’ 
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Fig. 2 The plots of approximate solutions of Example 2 for M=2. 


and the solutions of (l28l) can be found via (EB). The graphs of these solutions 
plotted in Fig. [B 


Example 3. Consider the problem of a half plane containing a crack par¬ 
allel to the boundary which illustrated in Fig. [3] and formulated as the system 

0 


[Kn(t,T)t/)i(T) + Ki 2 (t,T)(/) 2 (T)dT] = 0 

(33) 

^7^ dr + [K2l(t,T)(/)i(T) + K22(t,T)(f>2(T)dT] = TT 


with 


Kn(t,T) = 


K22{t,T) = 


T — t 

(r — -I- 4/1^ 


T — t 

(t — + 4:h'^ 


8h^(r — t) 

[(r - t)2 + 4 /i2]2 “ 
Ki2{t,T) = K2l(t, 
8h^(T — t) 

[(T-t)2+4h2]2 “ 


4/i2(t — t) \l2h? — (t — t)2] 
[(r - t)2 4h2]3 

8h^ [4^2 - 3(t - t)2] 
[(t - t)2 + 4h2]3 
4h^(T — t) [12/i2 — (t — t)2] 
[(r-t)2+4h2]3 


where h is the distance of crack from the boundary. The physical conditions 
of the problem impose that the relations 


(^i(T)dT = 0, 


02 ( t ) dr = 0 




( 34 ) 


and 

01 (t) = 01 (-t) , 02 (t) = -02 (-i) 


to be satisfied. Therefore the unknown functions may be expressed as 
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Fig. 3 Crack parallel to a free boundary 


M M 

~ ^ ^2j(r), (t>2{T) /i_^2 XI T2j-l{T). (35) 

V j=o ^ i=i 

For = 1, it follows from the orthogonality condition (USD that the second 
condition in (IMl) satisfies and the first one gives /?io = 0. 

By taking /3ii = ^21 = 0, as arbitrary values, the remaining coefficients jiij 
are uniquely determined from the linear algebraic system (HU) for each values 
of h and M. This leads us to find the functions 4’i{t) and from (( 551 ) . 
The stress intensity factors 


fci = lim \/l — 02 (t) 
k2 = lim \/l — 01 (i”) 

and their absolute estimation errors(Est.Err.) reported in table [T] Eor h = 00 
and Kij(t,T) = 0, from Eqs. (1^ and (fT^ - (|T^ the exact solutions of (1^ are 
obtained as 

^i(t) = 0, (j)2{T) = ;y=^- 

which give fci = 1 and ^2 = 0. This is shown in the last row of Table [TJ The 
table shows the rapid convergence of the results even for relatively small values 
of M. 


4 Conclusions 

We described a new idea of using Chebyshev polynomials for the numerical 
solution of system O- As applications of this idea, we have solved the simple 
examples 1 and 2 as a coupled system of singular integral equations of the 
first kind. In example 3, we studied a crack problem in solid mechanics and we 
reported the numerical results (see Table [Ij to show the efficiency and rapid 
convergence of the proposed method for all these kinds of problems. 
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Table 1 Stress intensity factors for the crack parallel to the boundary 


h 

M 

ki 

Est.Err. k. 

k2 

Est.Err. k 2 

0.2 

6 

4.878800637605022 

6.3e-14 

1.750099102171126 

6.2e-14 


7 

4.788277537335018 

l.le-14 

1.727809740547429 

4.1e-15 


8 

4.760729834685963 

4.8e-14 

1.719782910590219 

l.le-14 

0.4 

3 

2.607272141646415 

4.3e-15 

0.7745787927510580 

1.6e-16 


4 

2.594500911475041 

7.3e-15 

0.7266641783709941 

5.0e-15 


6 

2.594423234973139 

4.2e-14 

0.7376171346942053 

3.3e-16 

0.6 

2 

1.834057544899021 

l.le-15 

0.5664257041432605 

4.1e-16 


5 

1.960455689663461 

6.1e-15 

0.4297949760368867 

1.6e-15 

0.8 

2 

1.608371955353828 

1.6e-15 

0.3323260582700188 

2.8e-16 


3 

1.660617572058080 

8.7e-16 

0.2675691556476836 

4.6e-16 

1.0 

2 

1.461157081431933 

2.0e-15 

0.2104682299562445 

l.le-16 


4 

1.485914720666516 

2.1e-16 

0.1796691052492212 

l.le-16 

1.2 

4 

1.372176156193755 

5.0e-16 

0.1234414146531335 

0. 

1.5 

4 

1.262800608570183 

1.6e-15 

0.07465158121522054 

1.7e-16 

2.0 

3 

1.162112249974693 

l.le-15 

0.03662808437088003 

0. 

3.0 

2 

1.077621553329114 

3.3e-16 

0.01274529646673066 

0. 

10 

2 

1.007451045420713 

2.6e-16 

0.00037197952964307 

0. 

oo 

1 

1 

0 

0 

0 
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